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Abstract 

For 71=1 — exp(— 7?.) a random closed set obtained by exponential transformation of the closed range 
7?. of a subordinator, a regenerative composition of generic positive integer n is defined by recording the 
sizes of clusters of n uniform random points as they are separated by the points of TZ. We focus on 
the number of parts Kn of the composition when TZ is derived from a gamma subordinator. We prove 
logarithmic asymptotics of the moments and central limit theorems for Kn and other functionals of the 
composition such as the number of singletons, doubletons, etc. This study complements our previous 
work on asymptotics of these functionals when the tail of the Levy measure is regularly varying at 0+. 
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1 Introduction 



For a drift-free subordinator {St,t > 0) define TZ C [0, 1] to be the range of the multiplicative subordinator 
St = I — exp{—St),t > 0, obtained by the exponential transform. The gaps in TZ are the interval 
components of the open set [0, 1] \ TZ. The following construction of random compositions of integers has 
been studied in [3 El El E| • For each n, let -ui, . . . ,u„ be an independent sample from the uniform 
distribution on [0, 1], also independent of TZ. Define an ordered partition of the sample into nonempty 
blocks by assigning two sample points Ui and Uj to the same block if and only if the points fall in the 
same gap in TZ. A composition C„ of integer n is defined by the sequence of sizes of blocks, ordered from 
left to right. The sequence of compositions (C„) is a regenerative composition structure in the sense of 
[^. The regeneration property means that, conditionally given the first part of C„ is r, the remaining 
composition of n — r has the same distribution as Cn-r- 

The purpose of this paper is to prove central limit theorems and determine the asymptotic behaviour 
of moments for the number of parts Kn and some other functionals of the composition C„. We focus on 
the gamma subordinator with parameter 9 > 0, that is with Levy measure I'idy) = (e~^^/?/) dy ,. More 
generally, we obtain similar results for Levy measures which are sufficiently like the gamma Levy measure 
in their asymptotic behaviour at both 0+ and oo^. In the case of gamma subordinators our principal 
result fThcorcm lll|l specialises as follows: 

Theorem 1 For the gamma subordinator with parameter 9 > the number of parts Kn is asymptoti- 
cally normal, with Eif„ ~ (9/2) log^ n and vav Kn ~ (^/3) log"^ n. 

Other variables under consideration are Kn{t), the number of parts of a partial composition produced by 
the subordinator restricted to [0,t], and the counts Kn,r and Kn,r{t) defined as the multiplicity of part 
r in Cn and the multiplicity of part r in the partial composition, respectively. 

In our previous paper '10' we studied Kn for a subordinator whose Levy measure is regularly varying 
at with index a £ ]0, 1]. In that case we obtained limit theorems for Kn with normalisation by n" £(n) 
where ^ is a function of slow variation at oo. According to these results, _ftr„/(n" i{n)) converges strongly 
and with all moments to a nondegenerate limit, which is not gaussian. Moreover, Kn{t) is of the same 
order of magnitude as Kn, for each t > 0. For the case of a compound Poisson process, when the Levy 
measure is finite, Gnedin [Jj had previously shown that Kn typically exhibits a normal limit with both 
expectation and variance of the order logTi, though Kn{t) remains bounded as n grows, for each i > 0. 
Thus the case of gamma-type subordinators may be seen as intermediate between the above two: as in the 
regular variation case, Kn{t) is unbounded, but, as in the case of a finite Levy measure, its contribution 
to Kn is asymptotically negligible, for each t > 0. 

Results akin to Theorem ^ are very different from the limit theorems in IIUI and require other 
techniques. Our approach here relies on a recent version of the contraction method due to Neininger 
and Riischendorf |12j . Application of this method requires an appropriate decomposition of Kn, control 
of the principal asymptotics of two moments, and an estimate of the remainder term in the asymptotic 
expansion of the variance. We provide this background by poissonising Cn and applying the Mellin 
transform technique to analyse integral recursions for the moments. 

We also establish joint convergence to a gaussian limit for the sequence of small-part counts {Kn,r, r > 
1), as n oo. This kind of convergence holds neither in the regular variation case (when the scaled 
Kn,r^s converge to multiples of the same variable) nor in the compound Poisson case (when {Kn,r, n > 1) 
is bounded uniformly in n, for each fixed r). It resembles, however, Karlin's central limit theorem 
for nonrandom frequencies in the case of regular variation with a positive index. 

2 Setup and notation 

Each part of C„ corresponds to a nonempty cluster of the uniform sample points within some gap. So Kn 
is the number of gaps occupied by at least one of n sample points, and Kn^r is the number of gaps that 
contain exactly r sample points. Keep in mind that Kn = "Sj-Kn^r, and n = Tij-rKn^r- Sometimes Kn 
is called the length, and n the weight of the composition. Similarly, Kn(t) and Kn,r{t) are the counts of 
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sizes of clusters within the gaps of a smaher set TZ n [0, St], corresponding to the subordinator restricted 
to [0,t]. 

For ly the Levy measure of (St) let ly be the measure on [0, 1] obtained from v by the exponential 
transform y i— > 1 — e^^ . The Laplace exponent of the subordinator is given by the formulas 

/•oo /•! 

$(s)=/ (1 - e-"^)i/(d?/) = / {1 - {1 - xy)V{dx) , ^s>0. (1) 



^0 "'0 

Considering also the binomial moments 

^{n:m):^J^ (^^^ a;"(l - 2;)""'";7(dx) , 1 < m < n , (2) 

the ratio $(71 : r7i)/$(n) ,m — 1, . . . ,n, gives the distribution of the first part of C„, and the probability 
of a particular value of C„ is a product of factors of this type We also record the relation between 
power and logarithmic moments of measures and 1^: 

z.(dy) = / I log(l - x) P Didx) , j = 1, 2, . . . . (3) 

JQ 

As in |lUj . we shall also consider the poissonised composition Cp derived from TZ by the same con- 
struction as Cm but with the set of atoms of a homogeneous Poisson point process on [0, 1] with rate p, 
instead of the uniform sample of fixed size. We denote Kp, Kp^r, Kpit), Kp^r{t) the obvious counterparts 
of the fixed-n variables, and introduce the Poisson transform of the Laplace exponent 

1 °° 



/ (1 - e-n^dx) - ^ - <i>(n). (4) 

We shall be using throughout the abbreviations 

L = logp, log"a; = (logx)". 

Our primary concern are the compositions induced by a gamma subordinator, whose Levy measure, 
its exponential transform and the Laplace exponent are given in terms of a parameter 9 > hy 

Kdy) = ^dy, y e]0,oo[, ^dx) ^ ~f ' dx , xe]0,l]. - log (l + ^) , (5) 

y |log(l-a;)| V 9 J 

The density of ? has a pole and the tail Dlx, 1] has a logarithmic singularity at 0. We could consider a 
larger family of measures which differ from ((SJ by a positive factor, but this would not increase generality, 
since multiplying the Levy measure by such a factor does not affect the laws of TZ, Cm Kn or Kp (although 
it requires a linear time-change for functionals of a partial composition like Kn{t) or Kp{t)). 

More generally, we consider a drift-free subordinator with Levy measure ly, which has a continuous 
density on ]0, oo[ and satisfies the following conditions (L) and (R) which we record in terms of v and D: 

(L) either of the following four equivalent conditions holds 

(f>(p) = L + c + 0(p"^) , as p — > cx) , $(p) L + c + 0(p"') , as p — > 00 

iy[y, 00] = - log y + c - 7 + ©(y") , as y i , v[x, 1] = - log a; + c - 7 + 0{x'') , as a; i 
where c and e > are some constants and 7 — — r'(l) is the Euler constant, 
(R) either of the following two equivalent conditions holds 

iy[y, 00] = Oie^^y) , y T 00 , ^x, 1] = 0((1 - xY) , xU 

where e > 0. 
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See Appendix for the equivalence of conditions in (L) (L stands for left and logarithmic at 0) . Condition 
(R) (R for right or regular at 1) implies that all moments mj are finite, and that $ is analytical for 
5Rs > -e. 

Throughout in this paper, e or ^ are some sufficiently small positive constants whose values are 
context-related and may change from line to line. We denote c the constant in (L). For a generic positive 
constant we write d, and use Cj, dj for further real constants which are not important, but may appear 
with explicit evaluation in asymptotic expansions or other formulas. For shorthand, the right tail of v is 
denoted by 

v{x) — v[x, 1] , X <E ]0, 1] 
and a homogeneous Poisson point process with intensity p on [0, 1] is denoted PPP(p). 

3 Recursions 

We shall make use of two types of recursions. These apply to regenerative compositions generated by a 
drift-free subordinator, without any special assumptions on the Levy measure. (They can also be readily 
generalised to the case with drift). The first type of recursion, suited to the Poisson framework, is based 
on a decomposition of TZ analogous to the first-jump decomposition of a renewal process. The second 
type is based on splitting the range of a subordinator by its value at a stopping time. 

Integral recursions Define a pattern i? to be a nonempty set of positive integers. We say that a 
finite configuration of points within a given subinterval ]a,6[G]0, 1[ fits in E if the cardinality of the 
configuration is in E. For E fixed, let 

be the probability that the configuration of atoms of PPP(p) fits in the pattern E. The probability that 
the PPP(p) configuration on a subinterval ]a, b[ fits in E is equal to 7r(p(6 — a)). 

Consider the poissonised composition induced by an arbitrary drift-free subordinator with v{l} = 
(no killing). For E fixed, let Np be the number of gaps of TZ such that the PPP(p) configuration within 
the gaps fits in E. The count Np is a functional of the poissonised composition Cp\ in particular, defining 
E to be {1,2, . . .} (the configuration is nonempty) or E = {r} (the configuration consists of r points) 
we obtain Np — Kp, respectively Np — Kp,r- Other possibilities may be considered, for example taking 
E — {1, 3, 5, . . .} the count Np becomes the number of odd parts of the composition. 

Let Pj{p) — P(iVp = j) be the distribution of Np for some fixed pattern. Each pj may be extended to 
an entire function of a complex variable, with the initial value Pj{0) ^ l{j = 0). Introduce the factorial 
moments 

/(™)(p) - ENp{Np - 1) . . . (iVp - m + 1) , m = 0, 1, . . . 
with the convention f^'^\p) — 1. The following lemma is a minor variation of |l()l Lemma 6.1]. 

Lemma 2 Let E he a pattern with the probability of occurrence 7r{p). The distribution of Np satisfies 
the integral recursion 

(Pi(/') - (1 - 7r(pa;))pj(p(l - x)) ?(da;) = / 7r{px)pj^i{p{l - x))v{dx) (6) 

Jo 

for j = 1,2,... and p > 0. The same equation is also valid for j — 0, with on the right-hand side. 
Furthermore, the factorial moments of Np satisfy the recursion 

l\f^'^\p)-f^"'\p{l~x))V{dx)^m [\{px)f^^-'\p{l-x))D{Ax) (7) 
Jo Jo 

which taken together with f'^^^p) — 1 and f^™\0) = l{m — 0) uniquely determines them. 
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Proof. The derivation of recursions follows as in Lemma 6.1]. The uniqueness claim is a consequence 
of analyticity. □ 



Integral recursions for the distribution and the factorial moments of Kp follow by taking 7r(p) = 1 — e 
which is the probability that the PPP(p) configuration is nonempty; we have then po(p) — e^''- To obtain 
recursions for K^^p we should take 7r(p) = e^^p'' /r\ (in this case no simple formula for po{p) is known). 

Splitting at an independent exponential time Further recursions follow by splitting the range of 
a subordinator by its value at a stopping time. Though the fixed- rt version is needed to apply we 
focus on the poissonised model, which simplifies moment computations. Transfer of results to the fixed 
n model then follows by elementary depoissonisation. 

For each t >0, consider the sigma-algebra generated by both {Su,u < t) and the restriction of PPP(p) 
to [0, St]. As t varies, this defines a filtration, and we can consider stopping times with respect to it. Let 
T be such a stopping time with range [0, oo\ and let hp = Kpij) be the number of parts of Cp produced by 
the subordinator up to time r. The strong Markov property of TZ along with the independence property 
of PPP imply that the number of parts satisfies the distributional equation 

kp^hp + Kl (8) 

where bp = Kp{T) , Ip — p{l — Sr) and (K'p) is a distributional copy of {Kp), independent of {bp, Ip). For 
example, letting r be the first time that a jump of the subordinator covers at least one Poisson point, 
the equation holds with bp = 1. Identities analogous to ^ can be written also for Kp^r and more general 
pattern counts. 

We shall consider in some detail the most important case when r is an exponential time with rate 
A, independent of the subordinator and the PPP (thus r is a randomised stopping time with respect 
to the above filtration). The stopped process has an obvious interpretation as a killed multiplicative 
subordinator, which jumps at time r to the terminal value 1 (thus producing the final meander gap in 
the range). In accord with our previous notation, let Kp{T) be the number of parts produced by the 
subordinator before r (thus a possible part induced by the meander gap is not taken into account). Let 
{pj{p),j > 0) and {f^"^\m > 0) be the distribution and factorial moments oi Kp{T). 

Lemma 3 Forr an exponential time with rate \, independent of the process {K p{t),t > 0), the factorial 
moments of Kp{T) satisfy the recursion 

A/^™Hp)+ /\/(™)(p)-/(™Hp(l-x)))i7(dx) [\{px)f^^''Hp{l^x))^dx) (9) 

Jo Jo 

which taken together with f'^^^p) — 1 and f^™\p) — l{m — 0) uniquely determines them. 
Proof. In the case of a finite Levy measure, the first-jump decomposition of the range yields 

Pj^P^^ (^(1 - 7r(pa;))pj (p(l - x)) + 7r(px)pj_i(p(l - x)) 

, . A T'{dx) 

To see the extension for an arbitrary Levy measure, substitute 

A /"^ / N i^{dx) 

^^^p^ ^ p^^p^mTx + Jo ""^^'^mTx ' 

then rearrange terms and argue as in |1U', Lemma 6.1]. □ 



^{dx) 
i7(0) + A 
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If the meander gap is counted (this corresponds to Kp for killed subordinator) we should modify the 
recursion for moments by adding A7r(p)l(m = 1) to the right-hand side. 

Remark. The distribution and moments of Kp{t) may be obtained as the inverse Laplace transform in 
A from the analogous quantities for Kp{T). 



4 Moments of Kp 

Mellin transform resolution The recursion ^ is intrinsically related to a convolution- type integral 
equation 

ifip) - fipil - xWidx) = gip) (10) 

with the function g given and the function / unknown. For g analytic, there is a unique analytic solution 
/ with given initial value /(O). Observe that constants d constitute the null space of the integral operator, 
thus for each solution / the function / + d is another solution, with the initial value /(O) + d. 

It is easy to write out a power series solution to (IIUII , but such a representation itself does not help 
describe the large-p behaviour. We turn therefore to asymptotic methods based on the Mellin transform. 
Recall that for a locally integrable real- valued function (f) on [0, oo] the Mellin transform is defined by the 
integral 

/•oo 

M0(s) = / p'-^Mdp 



Jo 

which is assumed to converge absolutely for s in some open interval on the real axis, hence also for s in 
the open vertical strip based on the interval. The left convergence abscissa of M cj) is determined by the 
behaviour of (j) near 0, while the right convergence abscissa is determined by the behaviour of at oo. 
The analysis to follow is based on the formula 

M/(.) = ^ (11) 

which is valid in the common domain of definition of all ingredients. The formula follows by Fubini and 
a change of variable from 

p^-Mp / {f{p) - /(p(l - x))V{dx) = Mf{s) / (1 - (1 - x)-')v{dx) = Mf{sM~s) 
Jo Jo 

according to formula Q. Another important tool is the following correspondence between asymptotic 
expansion of a function and singularities of its Mellin transform. See for a fuller exposition of this 
technique. 

Lemma 4 Section 2] Suppose the Mellin transform M f of a function f is analytic in a strip 
a < 5Rs < b. If M f can be extended meromorphically through the right convergence abscissa in a larger 
trip a < 5Rs < b + e and has finitely many poles there, then each pole z and each term d{s ~ z)~'°~^ in 
the Laurent expansion of M f at z contributes the term 

to the asymptotic expansion of f at oo. The remainder term of the expansion of f is then 0{p^^^'^). 
Conversely, the asymptotic expansion of f at oo with such a remainder implies the termwise singular 
expansion of M f in the strip a < 5Rs < b + e provided that for some d > 

\Mfis)\^Oi\s\-'-') 

as |s| — > cxD in a strip 6' < 5ft s < 6 + e for some b' < b. 
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For the expectation f^^\p) := Eii'p the formulas (0) and Hll|l specialise as 

M/W(.) = ^ (12) 

where ^{s) is defined by Using Fubini and integration by parts we transform the numerator as 

M$(s) = / p'^^dp pe^P''iy{x)dx = v{x)dx I p'e'P''dp = 
Jo Jo Jo Jo 

T{s + 1) [ x-'-^i^{x)dx^~r{s) f x-'T^idx). 
Jo Jo 

We can now re- write formula (|12|) as 

\^(\-x)-''~^v{x)dx 

where we used (0) and ^ in the form 

$(-s : -s) = -s / x^""^^ v{x)dx , $(-s) = -s / (\ - xy^^ v{x) 
Jo Jo 



dx. 



Expectation Assuming (L) we conclude from that Mf^^^ is defined in the strip —1 < 5Rs < 0. 
We focus therefore on the meromorphic continuation of M/^^^ through the right convergence abscissa 
5Rs = 0. The same formula says that the poles of A//^^^ might be caused either by poles of F or by poles 
of $(— s : —s)/s, or by zeros of s)/s. We shall consider the three ingredients separately. 
For 3?s > — 1 the gamma function has a unique pole at 0, with Laurent expansion 

r(s) -7 + rfis + 0(s2) 

where di = 7^/2 + 7r^/12. With reference to Stirling's approximation, when u is a bounded real number 
and is a large real number 

\T{u + ±v)\r^ (27r)i/2|^;p-i/2e-^l''l/2 (14) 

uniformly in u. 

By assumption (R), the function 

-$(-s)/s= / {1 - x)-'-^ i7{x) dx 
Jo 

is analytic for Sfts < e, and its behaviour at complex infinity is regulated by the next lemma. 

Lemma ^ If v has a continuous density on ]0,cxd[ and satisfies (L) and (R) then 

$(s) ~ 7 log \s\ , as \s\ — > oo 
uniformly in any strip — e < 3?s < d, for e sufficiently small. 

Proof. The measure y v{dy) is finite, with exponential decay at infinity, thus for s = u -\- Iv the partial 
derivative 



— $(u + ±v) = 
du Jo 



'^yv{dy) 



is bounded for u > — e, uniformly in v. Hence + it;) — $(iu)| is uniformly bounded both in v and 
in u G [— e, d]. Thus it is sufficient to show the asymptotics for u = 0. But for the Fourier integral 

/"OO 

<i>(iw) — iv / e^"'^j^[y, oo] dy 
Jo 
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the desired asymptotics follows from the expansion (L) and smoothness of • , oo] by application of 01 
Section 3.1, Theorem 1.11] (with a further reference to T). □ 



By one further elementary lemma (see Appendix), condition (L) implies that $(— s)/s has no zeroes for 
5Rs < 0, hence by Lemma |5l there are no zeroes in the strip — 1 < SRs < e. The Taylor expansion at 
involves the logarithmic moments Q 

= nil + m2:J + + 0{s^). (15) 



s -2! "3! 

The integral 

x^'^^^iy{x) dx 



-$(-s : -s)/s = / : 
Jo 



converges for 3? s < and by assumption (L) this function is meromorphic for 3? s < e and in this 
half-plane has the unique pole at s = 0, where the Laurent series starts with 

— $( — S ; — s) 9 / X 1 , ^ 

^ = s"^ - (c - 7)s-i +d2 + o(l) 



where c?2 is some constant. The meromorphic extension is obtained from this formula, provided we 
substitute a suitable analytic function for the o(l) term. This follows by computing 

1 

(I loga;| + c - 7 + h{x))x^''^^dx 

with h{x) = 0{x^) and 3? s < 0, and noting that the integral of h is analytic for 9? s < e and it is bounded 
in this domain due to a maximum ridge on the real axis. 

It follows that M f^^^ is meromorphic in the strip —1 < 3fJs < e with a unique singularity at 0, where 
it has a triple pole. Putting the ingredients together and with a minor assistance of Mathematica we 
obtain the Laurent expansion 

Mf^{s) = -1.-3 +(^ + ^\ + d,s-' + 0(1) 
mi \mi 2mj y 

where 

-di — ^2 — C7 + 7^ cm2 TB.2 m3 



d.= 



mi 2mi 4mf Bmi^^ 

is a constant whose explicit value will not be used below. By Lemma 0] as p — > cx), the asymptotic 
expansion of the expectation is 

2mi \ 2ml ™i 



(p) = —L' + K^ + _ L-d3 + O(p-). (16) 



The decay condition on M f'^^^ required in Lemma^is satisfied, because this Mf^^^^ has exponential decay 
as \s\ — > oo in a strip about the imaginary axis due to 114() . Lemma^and because |$(s : s)\ = 0{\s\~^). 

Evaluating the right-hand side of (jT))- We will use Hll(l once again, to compute f^^\ The right-hand 
side of IZj) is evaluated with the help of the next lemma. 

Lemma 6 Assume (L) and (R) and suppose f is an increasing positive function of p £ [0, oo], which 
admits an asymptotic expansion at oo as a polynomial in L, with a remainder 0{p^'^). Then a similar 
expansion holds also for the function 

h{p)^ f\l~e'nf{p{^'^))Hdx) 
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and this expansion is obtained from that of f by replacing each L*^ with 



j=i Vj/ 

where c is as in (L). 
Proof. By the assumption 

k 

/(p) = 5]d,L^+0(p-^). (17) 

j=o 

Start with the case f — iJ^ . The integral is computed by the binomial expansion of (L + log(l — x))^ 
and termwise integration. The leading term does not depend on x and, in view of (L), it contributes 
L*^ (L + c) + 0{p~'^). The lower order terms are integrated using ||2J), with a remainder estimated by a 
constant multiple of 

L'' [ e-''^|log(l-x)|'=i;(dx), k>l. 
Jo 

The last integral is evaluated as 0{p^~'^) using integration by parts, using (L) and applying a standard 
Tauberian theorem on Laplace transforms. 

By linearity we obtain asymptotic expansion for the model integral I{p) which corresponds to the 
polynomial part of /, as in H17|l but with zero remainder term. One checks easily using (R) that for 
5 := the contribution to /(p) of the integral over [1 — 5, 1] is estimated as 0{p^'^). 

For a general / as in H17|l we again split the integral at \ — 5. For a; e [0, 1 — (5] we can apply expansion 
(I17|l to f{p{l — x)) with a remainder 0{p^'^), thus by (L) we have 

' V - e-nfipi^ - ^mdx) = Hp) + o{p-^). 

It remains to show that the integral over [1 — 6, I] is 0{p^'^) but this is easy because the assumed 
monotonicity yields the bound 

' (1 - e-nf{p{^ - ^mdx) < f{p5) i7(l - S) 
lis 

and this has the desired order by H17|) . assumption (R) and our choice of 5. □ 



Taking f^^'^ for /, the lemma translates the expansion into the p ^ oo formula for the right-hand 
side of lO with m — 2 

gi{p) = -L' + + d,L + d, + 0{p-') (18) 

mi \mi m^J 

where dn, ds are some constants, and the factor 2 in lO is included in the definition of 171. 

Moments of the second order Using the direct correspondence in Lemma 0] we see that is the sole 
singularity of Mgi, and derive the expansion at s = 

M,i(.).l.--f^ + ^).-3^0(.-) 
mi Vmi mf / 



which together with (|11() and H15() implies 



<P(— s) mf mi 
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Finally, the inverse correspondence applied to Z*^^-* yields 



= -l^L^ +(^ + ^)L^ + dL'+ 0(L) (19) 



4m2 \ml 6mlJ 

provided the bound |M/(2)(s)| = 0{\s\-^-') can be established. This requires some effort because for 
Mgi no explicit formula is available. Postponing justification of the decay condition, and recalling 

EK, = f(^\p) , vari^p = f(^\p) + f^^Hp) - if^'Hp))' 
we have from (|16|l and p9|l 

Theorem 7 If the Levy measure of subordinator has a continuous density on ]0,oo[ and satisfies the 
assumptions (L) and (R), then as p ~* oo with L = logp, the asymptotic expansions of the first two 
central moments of Kp are 

IEi?p = -^l2 + 0(L) (20) 
z mi 

var K =p^lJ^ + 0(l2) . (21) 



Justification of the decay condition Our plan is to decompose gi = /ii + h2 so that hi will absorb 
the logarithmic part of gi at oo and will have a manageable Mellin transform, while M/i2 will satisfy the 
decay condition by virtue of the following classical result. 

Lemma 8 jl6l Section 1.29] Let (j) be analytic in a sector —a < Argp < /3 with < a, P < tt, and 
suppose that in this sector for some a < b and each 5 > 

(^{p) = Oi\pr-') as IpUO, 0(p)= 0(1^1-"+^) as |p| T oo . 

Then M(j} is analytic in the strip a < ^ s < b and satisfies 

M0(s) = O (e^^'^-'^)'*") as5s^cx3, M0(s) = O (e^""'^)'*") as 9s ^ -oo , 

for each e > uniformly in every strip strictly inside a < 5Rs < 6. 

By definition 

[' f^'Hp{l-x)Kl-e-n^dx) 
Jo 

which is an entire function, where 

/W(p) = e-''f;^EX„ 

n=l 

and the series in p has all coefficients positive. The same applies to the series involved in 
Therefore \gi{p)\ < gi{\p\) and by ^ 

|5i(p)|=0(L3) ^oo. 

Writing as 

3 

gi(p) = ^djL^+0(p-^), p^oo,peM 

3=0 
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and selecting 

3 

h^{p) := doe-'/p + log(p + 1) ^ d, log^-i p 

we obtain a function h2 '■= gi — hi which is analytic in the open halfplanc p > 0, has the same expansion 
in powers of L as gi for p ~^ oo and satisfies Lemma |S1 with a = I, b = e and a = f3 = tt/2. Hence by the 
lemma 

\Mh2{s)\ = O (e-(''/2-5)|s|\ ^ fQ^. _ 1/2 < gfjs < e , |s| ^ oo 



for each d. The meromorphic continuation of Mhi from — l<5Rs<0to the halfplane 5Rs > — 1 follows 
from the standard Mellin transform formulas 

Me-'/ ■ (s) - r(-s), A/log(l + ■){s) = -T{s)T{-s) , M(0( • ) log( • = 

which by application of (|14() also imply 

|Af/ti(s)| = O (e-('^/2-*)l''l) , for - 1/2 < SRs < e , |s| ^ c5o . 



It follows that in a strip containing the imaginary axis |Mgi(s)| = O (e '^'^^^ ^)\'^\^^ and the same 
estimate holds for |M/(^)(s)| = |A/gi(s)/$(— s)| due to Lcmma|Sl Thus \Mf^^^\ decays exponentially 
fast and Lemma 0] is applicable. 

Number of parts prior to the exponential split Let r be an independent exponential time, with 
rate A. By the Mellin transform in p satisfies 

M{EK,{t)){s)- 



A + $(-s) 



The term A in the denominator kills the zero of <& at s = 0, therefore arguing as in the analysis of (|13|l 
we see that this expression has only double pole at s = 0. It follows that 

EKp{T) =dL + 0(l) 

for each A > 0, with d depending on A. Repeatedly appealing to Q and Lemmas HJ and [S] we obtain 
j(m)(p) ^ 0(L'"), which trivially impHes the following rough estimate 

Lemma 9 For t an independent exponential time, as p oo with L = log p, 

E(i?p(r))'" = 0(L'"), m=l,2,... 



Depoissonisation Modifying the argument in |l()l Section 6.2] we obtain for p = n ^ oo 

for each i > 0, almost surely, with all factorial moments, and with at least two central moments. The same 
sandwich-type proof works because the moments grow logarithmically (thus vary regularly as n ^ oo). 

5 Central limit theorem for Kn 

We switch to fixed-n framework, in order to apply a delicate recent result due to Neininger and Riischendorf, 
which we reproduce below with a minor adaptation and notational changes. Suppose a sequence (Yn) of 
random variables satisfies 

Yn = bn + 17„ , n>l (22) 

where {Yn) = iX^^), the pair {bn,In) is independent of (1^), each /„ assumes values in {0, . . . and 
P(/„ = rt) < 1 for n > 1. Let p„ = EF„, cr^ — vary„, and || • II3 denote the £'^-norm of a random 
variable. 
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Theorem 10 ^1 Theorem 2.1] Assume that each \\ Y„ ||3< oo and that {Yn) satisfies the recursion 
Suppose that for some constants C > 0, a > the following three conditions all hold: 



(i) 



lim sup E log 



/„ V 1 



< , sup 

n>l 



log 



/„ V 1 



< oo 



(ii) for some A G [0, 2a[ and some k 

II 6„ - M« + M/„ lis- Oaog" n), al^C log'" n + 0(log^ n) 

(ui) 



1 

a > — + max I k , 



A 



Then the law of 

VClog^n 

converges weakly to the standard normal distribution. 



Condition (i) says that the split In/n must be bounded away from and 1, and the further two conditions 
require logarithmic growth of moments. See [12) for proof and estimates of the distance between the 
probability laws. 



To apply the result to Kn we split the range of subordinator at an independent exponential time r 
with mean 1. Let 6„ = i^„(r) be the number of parts of C„ produced by the multiplicative subordinator 
in the period between to r. Observe that H22|l holds, with J„ being the number of uniform sample 
points larger than Sr ■ The asymptotics of moments in Theorem |3 and Lemma El apply to Kn and 5„ 
literally without change, by the virtue of depoissonisation, 

Let us check now the conditions of the theorem. Since 1 < Kn < n, all absolute moment^ of Kn 
are finite. Furthermore, the variable In/n converges strongly and with all moments to 1 — Sr, thus 
— log(/„/n) approaches the stopped value Sr of the additive subordinator, which is positive and has 
E,Sr < OO as a consequence of ma < oo, by an easy application of the Levy-Khintchine formula. Whence 

(i) is satisfied. In evaluating the asymptotics of moments we can switch to the poissonised composition 
with p — n, Ij — logn. Then from Theoremdwe see that || /i„ — Ha is of the order of L. By Lemma 
l^lalso II bn ||3= 0(L), because = 0{L^). By the above and Theorem[3the parameters involved in 

(ii) are a = 3/2, A = 2 and k = 1, so condition (iii) is satisfied. Thus all conditions of Theorem 1101 are 
fulfilled and we deduce the following conclusion: 



Theorem 11 If the levy measure has a continuous density on ]0, oo[ and satisfies (L) and (R) then 
the distribution of the random variable 

Kn - LV(2mi) ^ , 

L = log n 



v/ni2/(3m3) L3/2 ' 
converges weakly to the standard normal distribution, as n 
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6 Small parts of the composition 



We have shown in 10 that in the regular variation case all Kn,r, r — 1,2, . . . are of the same order of 
growth as Kn and, suitably normalised, converge almost surely to multiples of the same random variable. 
In the compound Poisson case these variables are bounded as n grows 7 , and if the increments of (St) 
are exponentially distributed (this corresponds to the Ewens partition structure) the Kn^s converge to 
independent Poisson random variables Our next goal is proving a joint central limit theorem for the 
small-part counts. 



Marginal central limit theorems The line of argument repeats that for Kn . Let /i™'' (p) be the mth 

factorial moment of Kp^r, e.g. f[^\p) is the mean number of singletons. Recall that the recursion lO 
holds for Z*^™^ — fr"^^ with 7r(p) = e^^p^ /r\. Applying the Mellin transform to the recursion we obtain 

Mfi^\s)^^±±^^hl^, -l<5i.<0. (23) 
r! <P(— s) 

This agrees with (|13f) and Kp = -^p,r in view of the identity 

r— 1 ^ ^ 

The right-hand side of ^ can be extended meromorphically through the imaginary axis, with a double 
pole at 0, where we have the Laurent expansion 

M/«(s) = — s-'-dis-i + 0(l) (24) 

where 

2cmi - 27 ml -I- m2 - 2mi (r) 
2nif r 

and '0^'^^ is the digamma function. The decay at complex infinity is justified as before, thus by Lemma^l 

f^'Hp) = — + di+0{p-^). (25) 

Now we need to translate (|24f) into asymptotics of the right-hand side of (TJ with m — 2. We note 
that evaluation of integrals with a factor of e~''^x^ essentially amounts to Tauberian-type asymptotics, 
since for r > 1 the factor compensates the singularity of at x = 0, and e"''^ is negligible outside 
each fixed vicinity of 0. In particular, we have 

Lemma 12 // the condition (L) holds then for p — > 00 



5(d., = i + 0(.-), 



r! ' r 



Proof. Integrating by parts and replacing v by the right-hand side of (L), and pushing the upper inte- 
gration limit to 00 the claim is reduced to the standard integral 

'e^P^p^x'^^^ e~'^^ p'^'^^x^ \ 1 
{- \ogx + d) [ n = I 







(r-1)! 



which does not depend on d, because 

p-px r r-l -px r+1 r\ 

e px p ^ 



(r-1)! 
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as is easily checked. 



□ 



Using the lemma and (|25|l we compute the right-hand side of Q as 

, , 2L 2di ^, , 
5P = — + — + 1 , 
mir^ r 



hence by Lemma 0] 



We compute then 



which yields again by Lemma ^ 



(the decay condition is again checked by application of Lemma Isj. This together with 1)2 5() implies 

varif,p= ( ^ + — )L + 0(1). 
\r^mj mi / 

Decomposing as in © at rate 1 exponential time, the Mellin transform of the expectation of hp — 
Kp^r(j) has <&(— s) + 1 in the denominator, hence all moments of hp remain bounded as p ^ oo. Passing 
to the fixed-71 version and applying Theorem II 01 with a — 1/2, A = k = we obtain 

Theorem 13 Under our assumptions on the Levy measure, for large n the distribution of each K^^r 
is approximately normal, with moments 

EKn,r = — logn +0(1), varfi:„,,. = ( ^ + — ) logn + 0(1) . 
mi?' \r^m^ rmi / 

The same is true for Kp^r as p — > oo. 



Joint central limit theorem Our strategy to prove a joint central limit theorem for the small parts 
counts is to consider more general finite patterns and functionals of composition such as E,. a^-Kp j.. By 
linearity, each variable of this kind decomposes as in ©, hence the method we applied to Kp ^ can be 
used to justify the normal limit. As an instance of such a functional consider 

Kp,[r] '■= KpS + • • • + Kp^r 

and define f^^\p) to be the mth factorial moment of Kp^^j.]- The corresponding pattern is {1, . . . ,r}, 
thus the moments satisfy ^ with 7r(p) = Y^^j=i ^ ''P'^/j- Obviously from H25() 

/ff = ^L + 0(l) 
' ' mi 

where hr := X]j=i ^^'^ the harmonic numbers. Letting Mellin's machine roll to produce f^^^ 9 

(2) (2) 

Mg — > ^'^f [r] ^ J[r] variance is computed as 

varif,H=f^ + ^)L + 0(l). 
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Similar computation for the pattern E — and Theorem 1 1 31 yield the covariance 

cov{Kp,,,Kp,,)^ (^l, + l{t=j) — ] L + 0(1). 

Theorem 14 Under our assumptions on v, as n ^ oo, the infinite random sequence 

{{Kn,r - Eif„,,) log-1/2 „ ^ r^l,2,...,n) 
converges in law to a multivariate gaussian sequence with the covariance matrix 

+ . (26) 

Proof. To justify the joint gaussian law it suffices to establish convergence to a gaussian limit for each 
finite linear combination 



From the above facts about Kn,r and from H26() it follows that both the expectation and the variance of 
Qn grow like L. Also, due to an obvious additivity, satisfies a distributional equation of the type (|22|l 
provided we split the range of subordinator by the value at an independent exponential time r. Then 6„ 
is equal to the contribution to Qn by the jumps of subordinator before r. The conditions of Theorem I lUI 
are checked exactly as for Kn or Kn,ri thus by this theorem Qn is approximately gaussian. □ 



7 The gamma case and further examples 

In the case of the gamma subordinators we have 

c = -log0, m, = ^^^, 

and the singular expansion 

4>(-s : -s) = ~s~^ + (c - 7) - d2S + O(s^) 

with constant 

r\ x{l^x)0-^+\og{l~x) ^ 

''^L — -.log(i-x) — 

which does not seem to simplify. The Mellin transform of EA'p is given by the formula 

\ogil-s/0)i -log(l-x) 

which defines a function which is meromorphic for 5R s > — 1 with a sole singularity at s = 0. The Laurent 
expansion of M/^^^ at s = involves 

(7-1)! 72 
m, = . ' , c = -\og9, di + — . 

Leaving only principal terms, the asymptotics of moments is 

6 9 
EKp = + 0(L) , var Kp^ - + 0(L2) 

as was promised in Theorem ^ 
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Further examples Another instance of a gamma-type subordinator was introduced in [H] to describe 
a composition resembling the ordered Ewens sampling formula. This subordinator has (multiplicative) 
Levy measure 

u{dx) = {1 - xf^^x^'^ dx , xe]0,l] 
with parameter 9 > 0. The Laplace exponent given by the formula 

' fr{\j + o-i 3 + e-i + sJ 

is a function interpolating the generalised harmonic numbers 

which can be seen as a combinatorial analogue of the logarithmical Laplace exponent for the gamma 
subordinator. 

The basic characteristics of subordinator are readily expressed in terms of the polygamma function 
Thus we have 

v{x) = -\ogx-'ip^°\0) + O{x~^) 

as is shown by expanding 



,7{x) = / (1 - zf'^z^^ dz = - \ogx + c{e) - 7 + ci{e)x + ... , 

J X 

substituting 6* = 1 to see that c(l) ~ — — V''°-'(l)i then differentiating in 9 and sending x to obtain 

c'{9)^ ( {I- zf-^z-^\og{l- z)dz^-il;^^\9), (28) 
Ja 

whence c{9) = -^l;<'°\9). It follows as in Appendix or directly from (|27|l that the expansion at infinity is 

$(p) -L-^(0)(^?)+O(p-i). 
The moments are computed by further differentiating H28|l as 

/ {l-xf-^x-^\\ogil-xWdx = {-iy+^4!'-^'>{9). 



See |2| for computation of some densities related to this family of subordinators, and | 15| for further 
examples of Levy measures with logarithmic singularity. 

Oscillatory asymptotics for a discrete measure The constant term c assumed in (L) cancels in the 
asymptotics of the moments. However, this assumption is essential by our approach. Let us assume a 
bounded oscillating term in place of c, and examine how the asymptotics could be affected. For ease of 
computation we consider the atomic measure 

oc 

^'(da;) = ^Si/^j{dx) . 

For this measure $(n) is equal to the expected maximum in a sample of n geometric random variables, 
which was analysed in IBl Example 12]. 
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Denoting [ • J , respectively { • }, the integer and the fractional parts of a positive number, we have 

i7{x) = [— log xj = — logx — {— logo;} , X e]0, 1] 

thus there is a logarithmic singularity but the expansion (L) does not hold, because the second term 
oscillates between and 1. Condition (R) is satisfied since (^{x) = for a; > 1/e and all moments are 
finite, 

oo 

m,=^|log(l-e-^")|\ 
e.g. mi = 0.6843, m2 ~ 0.2345 (truncated at four decimals). Furthermore, 

It is seen that M$ has a double pole at s = with Laurent expansion 

Af $(s) ^1- 2±1I1 + 1(1 + 67 + 672 + n^) + 0(s) 
s"^ s 2 

and infinitely many simple poles on the imaginary axis at 

Sfc = 27rifc, keZ\{0}. 

The conclusion of Lemma 0] still holds, which can be justified by applying the inverse Mellin transform 
formula and integrating over increasing rectangular contours, as fc — > cxi, with sides 3s = 27rfc + l/2, 3s = 
-27rfc - 1/2, SRs = 1/2, 5Rs = -1/2. Thus 

$(p) =L + 7 + l/2 + 0(L) + C)(p-i+^), p->oo (29) 
where the contribution of the imaginary poles amounts to the term 0(L) given by the formula 

cf>{u)^~ r(2^ifc)e2-i'=«, 

feGZ\{0} 

which is a periodic function with period one and a tiny amplitude. The fluctuations of the 0(1) term in 
are not asymptotically negligible, though they are very small due to the fast decay of the Fourier 
coefRcients, e.g. 

r(si) = r(s_i) = (0.126 + 0.501 i) 10""^. 

(truncated at seven decimals). 

It follows that Mf^-^'' has a triple pole at and simple imaginary poles Sk — 27rifc, hence 

where is another periodic function. Thus oscillation prevails only in the third term in the expansion 
of f^^\ However, computing gi, as in Lemma|51 we should include L?^{p), thus arriving at 

51 (p) = — L3 + 02(L)L2 + ... 
mi 

with oscillating 02- As a consequence we will have 

1 L4 + ^3(L)L3 + ..., 
4m 

and this suggests that the principal 0(L"^)-term in the asymptotic expansion of the variance var Kp will 
oscillate, though we could not establish this fact rigorously. 

In a similar situation of sampling from the geometric distribution, the expectation of the number of 
different values in a sample also involves an oscillating second term, and the same is true for the factorial 
moment of order 2 However, this has no impact on the principal asymptotics of the variance: the 
oscillating terms cancel and the variance converges to a constant, as had been shown long ago by Karlin, 
see m Example 6, p. 385]. 
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8 Appendix 

Lemma 15 The four conditions in (L) are equivalent. 

Proof. The equivalence of expansions of v and u follows by the change of variables x = 1 — and 
y = — log(l — a;). For example, assuming the expansion of u we obtain that of u by using 

l/lx, 1] = log(l — x),oo] 

and substituting 

- log(- log(l - x)) = - log(x(l +x/2 + ...))=- logx + 0{x). 
The expansions of integrals are equivalent to the expansions of measures by writing 

POO pi 

Hp)/p= / e-^V[y,oo]d?/, Hp)/p= / e-''^u[x,l]dx 
Jo Jo 

and using the classical integral 

POO 

p e-pyiog{l/y)dy = L + j 
Jo 

together with standard properties of the Laplace transform. □ 
Lemma 16 The function 

$(s) = [ (1- (l-x)")P(da;) 
Jo 

has no zeros for Res > 0. And if $(s) = for purely imaginary s, then v is atomic, with support 
{1 - a'',k ^ 1,2, . . .} for some < a < 1. 

Proof For Us > and x S ]0, 1[ we have |(1 - a;)^| < 1. Therefore 3? (1 - (1 - .x)^) > and the integral 
cannot be zero. For real r the equality $(ir) = is only possible when 1 = (1 — x)"" holds i^-almost 
everywhere, but such x is of the form x = 1 — exp(— 27rfc/|r|) for some k = 1,2, □ 
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